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Abstract. We present calculations of elliptic flow and azimuthal dependence of correlation radii in the 
ellipsoidally symmetric generalization of the Buda-Lund hydrodynamic model of hadron production in 
high-energy nuclear collisions. We compare them to data from RHIC by simultaneous fits to azimuthally 
integrated invariant spectra of pions, kaons and protons-antiprotons measured by PHENIX in Au-fAu 
reactions at center of mass energy of 200 AGeV. STAR data were used for azimuthally sensitive two-particle 
correlation function radii and for the transverse momentum dependence of the elliptic flow parameter V2. 
We have found that the transverse flow is faster in the reaction plane then out of plane, which results 
in a reaction zone that gets slightly more elongated in-plane than out of plane. The model parameters 
extracted from the fits are shown and discussed. 

PACS. 25.75.-q, 25.75.Gz, 25.75.Ld 



1 Introduction 

Important information about the properties of extremely 
hot strongly interacting matter comes from the observa- 
tion of azimuthal anisotropies in non-central ultra-relati- 
vistic nuclear collisions. The second order Fourier compo- 
nent of azimuthal hadron distributions is connected vi^ith 
the azimuthal dependence of transverse collective expan- 
sion velocity of the bulk matter [Tl[2] . That is in turn deter- 
mined by the differences of the initial pressure gradients in 
the two perpendicular transverse directions, as well as by 
the initial geometry, the initial velocity and temperature 
distributions of the fireball, and the equation of state [31 
|4] . The anisotropic shape of the fireball measured with the 
help of correlation femtoscopy [5] at the instant of final 
decoupling of hadrons bears information about the total 
lifespan of the hot matter: with time the originally out-of- 
reaction-plane shape becomes more and more round and 
may even become in-plane extended W. Unfortunately, 
in determining the elliptic flow and azimuthally sensitive 
correlation radii individually two effects — spatial and flow 
anisotropy — are entangled. For example, the same ellip- 
tic flow can be generated with varying flow anisotropy 
strength if the spatial anisotropy is adjusted appropri- 
ately [7]. 

In general, the precise way of the interplay between the 
two anisotropies is model dependent. It has been studied 



and shown to be different within the Buda-Lund model 
[H] as well as the Blast Wave model [7 . 

In this paper we analyze for the first time azimuthally 
sensitive Hanbury Brown - Twiss (HBT) radii, using data 
from non-central heavy ion collisions within framework of 
the Buda-Lund model. Note that the model successfully 
describes data from central Au+Au collisions at RHIC, as 
measured by BRAHMS, PHENIX, PHOBOS, and STAR 
collaborations, including identified particle spectra and 
transverse mass dependent HBT radii as well as the pseu- 
dorapidity distributions of charged particles. The model 
was shown before to describe the transverse mass and 
pseudorapidity dependence of elliptic flow of identifled 
particles at various energies and centralities in ref. |S] . The 
Buda-Lund model formalism for non-central collisions, in- 
cluding elliptic flow and azimuthal angle dependence of 
HBT radii has been proposed first in [S]. The model is 
defined with the help of its emission function. In order to 
take into account the effects of resonance decays it uses the 
core-halo model [TU] . In the present study, we improve on 
earlier versions of the Buda-Lund model, by scrutinizing 
the various components using azimuthally sensitive HBT 
data. Eventually we utilize a model that includes as a 
special case of T.S. Biro 's axially symmetric and acceler- 
ationless exact solution of relativistic hydrodynamics [TT] , 
in contrast to the original, earlier variant, ref. [5], which 
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was based on an ellipsoidally symmetric, but also non- 
accelerating exact solution of relativistic hydrodynamics, 
given by ref. [TJi ■ Similarly to ref . ;13, , we present an im- 
proved calculation, using the binary source formalism, to 
obtain the observables by using two saddle-points instead 
of only one. This results in an oscillating pre-factor in front 
of the Gaussian in the two-pion correlation function that 
we take into account for the formulae of the HBT radii. 

Azimuthally sensitive HBT radii were also considered 
recently in cascade models, e.g. in the fast Monte-Carlo 
model of ref. [TS], or, in the Hadronic Resonance Cas- 
cade [TB]. 

Data analysis of correlation HBT radii performed ear- 
lier with the Blast Wave model indicates that the fire- 
ball at the freeze-out is elongated slightly out of the re- 
action plane |17| . i.e. spatial deformation is similar as in 
the initial state given by the overlap function. This is also 
supported by the theoretical results from hydrodynamic 
simulations [6,18j and URQMD 19[. It sets limitations 
on the total lifespan. From all previous analyses it seems, 
however, that the final state anisotropy has an interest- 
ing non-monotonous dependence on collision energy with 
a minimum at the SPS energies [19J. In our analysis of 
the same data with a different model we observe for the 
first time at RHIC an in-plane elongation of the fireball 
at freeze-out. 

The paper is structured as follows: In Section [2] the 
basic features of the ellipsoidally symmetric Buda-Lund 
model are summarized. In Section|2]we derive the analytic 
formulae for the observables such as elliptic flow and the 
azimuthally asymmetric correlation radii. In Section 2] we 
show the results of the simultaneous model fits to exper- 
imental data from non-central collisions and we compare 
them to the ones obtained from fits to central data. In 
Section [5] our conclusions are presented. 



2 Buda-Lund model: basic features 

We restrict ourselves here to a short description of the 
model, for details see refs. [3120] • 

In the Buda-Lund model, the emission function is given 
by that of a hydrodynamically expanding fireball (core), 
surrounded by a halo of long lived resonances. The core 
emission function looks like: 



B{X,P) + Sq 



{I) 



where g is the degeneracy factor {g — \ for identified pseu- 
doscalar mesons, g = 2 for identified spin=l/2 baryons), 
and p^d'^E^lx) is a generalized Cooper-Frye term, de- 
scribing the fiux of particles through a distribution of 
layers of freeze-out hypersurfaces, B{x,p) is the (inverse) 
Boltzmann phase-space distribution, and the term Sq is 
determined by quantum statistics, Sq = 0, —1, and +1 
for Boltzmann, Bose-Einstein and Fermi-Dirac distribu- 
tions, respectively. Note that x^ = {t,rx,ry,rz) and p^ — 
{E,Px,Py,Pz) are the four- vectors of the space-time point 
x and the momentum p. 



For a relativistic, hydrodynamically expanding system, 
the (inverse) Boltzmann phase-space distribution is 



B{x,p) = exp 



p ■ u{x) ^{x) 
T{x) ~T{x) 



(2) 



The forms of the fiow four- velocity {ui,{x)), chemical po- 
tential {fi{x)), and temperature {T{x)) distributions are 
introduced below. Note that it can be mapped onto ex- 
act solutions of hydrodynamics, both in the relativistic 
and in the non- relativistic cases, as detailed in ref. |20) . 
For example, let us mention, that in the non-relativistic 
limit, the Buda-Lund hydro model corresponds to the ex- 
act, parametric, ellipsoidally symmetric solutions of non- 
relativistic hydrodynamics in ref. [4j which solution at late 
times converges to an accelerationless exact solution of rel- 
ativistic hydrodynamics as detailed in ref. I21| . According 
to our best knowledge, no similar connection has been ex- 
plored yet in case of the azimuthally sensitive version of 
the Blast Wave model of ref. , and known exact para- 
metric solutions of (relativistic) hydrodynamics. 

The generalized Cooper-Frye pre-factor, as described 
in refs. [^1^ was 



p'^d^S^ix) ^p^u^{x)H{T)d^x. 



(3) 



The time dependence of the emission, described by H{t) 
was approximated with a Gaussian distribution around 
the freeze-out proper-time tq. 



H{r) 



1 



(27rzir2)i/2 



exp 



(4) 



with At being the duration of the particle production 
in longitudinal proper-time (r = -^/t^ — r^). Of course, 
this function H(t) can be easily generalized to have more 
complicated forms, but the data discussed in the present 
paper do not require us to go beyond the Gaussian ap- 
proximation. However, we found that the analysis of the 
azimuthally sensitive HBT radii was actually sensitive to 
the structure of the Cooper-Frye pre-factor. Eq. ([3]) corre- 
sponds to freeze-out hypersurface layers that are pseudo- 
orthogonal to the four-velocity. For fiow profiles with sig- 
nificant longitudinal and radial fiows, as specified below, 
these hypersurfaces have positive correlations between the 
transverse radial coordinates rt and time t. When defin- 
ing the axially symmetric Buda-Lund model in ref. jl3) , 
such positive (rj , t) correlations were neglected and the 
freeze-out hypersurface was assumed to be a constant in 
the transverse direction. Recently, new exact analytic so- 
lutions of relativistic hydrodynamics also lead to freeze- 
out hypersurfaces with the property of nearly negligible 
{rt,t) correlations, see ref. Based on favourable com- 
parisons with data and analogies to the axially symmetric 
Buda-Lund model, we decided to keep this kind of freeze- 
out hypersurfaces for the purpose of the present paper. 
Note also that following ref. [13] of the axially symmetric 
case, we include a factor r to 7J*(t) and approximate it 
by a Gaussian H{t). Thus our modified Cooper-Frye term 
reads as 

p^d^Uf^{x) = rrit cosh(?7 — y)H{T)dT tq drj dr^dry, (5) 
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where rrit = \/ + -pi is the transverse raass, Pf is the 
transverse momentum, y is the rapidity and r\ is the lon- 
gitudinal space-time rapidity r\ — 0.51og[(t-|-rz)/(< — Tz)]. 
The four- velocity field, u'^{x) assumed to be a directional 
Hubble flow, where two different radial Hubble constants 
{H^ and Hy) characterize the different strength of the flow 
in the impact parameter plane {x—z) and out of this plane 
(y — z). This transverse flow is assumed to develop on a 
flow profile that is Bjorken-type at the rt = axes of the 
collision : 

u^{x) = (cosh[ry] cosh [jyt ], iJ^jr^;, iJj^r^^, sinh [77] cosh[?7t]) 

(6) 



exponentials in the pre-factor generate two terms with sep- 
arate saddle points. Hence, in the final formulae only x ap- 
pears. The saddle point coordinates — (r cosh(77),fa;,ry, 
rsinh(?7)) and the longitudinally boost invariant average 
emission widths (zir, Z^?;, Ry) are given as: 



To, 



yo-y 



1 + Arfmt/To 



PiRiRi/To 

lor I 



l + {a^ + Rl)E/To 



x,y, 



(9a) 
(9b) 

(9c) 



with Hx = Rx/Rx and Hy = Ry/Ry are the defining 
relationships for the flow in the impact parameter plane 
(called also as reaction plane) and in the remaining orthog- 
onal transverse direction. R^ and Ry are the characteristic 
geometrical system sizes in the two transverse directions, 
whereas Rx and Ry are their time derivatives. The aver- 
age transverse fluid rapidity ijt is also introduced with the 
defining relation sinh^[77t] = r^H^ + r^Hy, which ensures 
that u ■ u = 1. Note that at mid-rapidity, with 77 = 0, 
the velocity profile of T.S. Biro 's axially symmetric and 
accelerationless exact solution of relativistic hydrodynam- 
ics |TT] corresponds to the Rx/Rx = Ry/Ry case, while the 
solution discussed in ref. [12] coincides with this velocity 
profile. 

For the fugacity distribution exp we assume 



T{x) 



To 



2RI 



2RI 



so that it leads to a Gaussian in coordinate space. Arj 
denotes the space-time rapidity width and yo is the mid- 
rapidity. 

For the temperature profile we use the following form: 



1 



r(x) 




where Tq is the temperature of the collision center at the 



control the transversal and the temporal changes of the lo- 
cal temperature profile. Note that its dependence on the 
coordinate 77 will not be studied here (i.e. a,, — is as- 
sumed). 



3 Observables from the Buda-Lund model 

The observables can be calculated analytically from the 
Buda-Lund hydro model, using a double saddle-point ap- 
proximation in the integration. We quote results from ref. 
PU] . Note that in the binary source formalism of ref. [13] 
(Section 8 and 9) the double saddle-points are generated 
from the saddle-point x of the Boltzmann term by the 
product with the Cooper-Frey pre-factor in which the two 



At 



Ri — - 



1 + alE/TQ 

A,f 
" 1 + ArfE/To ' 

Ri 

l + ia^ + Rf)E/To 



) for j = X, y, 



(9d) 

(9e) 

(9f) 
(9g) 



where E = rrit cosh(^ — y). The invariant momentum dis- 
tribution is evaluated using an ellipsoidally symmetric gen- 
eralization of eqs. (127, 130-140) of ref. [13], that were first 
derived for the case of axially symmetric collisions: 



E- 



(7) where 



Ni{p) 



V 



:EV C 



1 



B{x,p) 



,0/0 At 

{2nf^ — R^^RxRy 



exp(Z\77^/2). 



(10) 

(11) 
(12) 



In the latter two expressions we use the notation = 
ToAfj and A* is the {y and pt dependent) intercept param- 
eter of the two-particle correlation function. 

The axially symmetric limit of these formulae, given in 
ref. [12 corresponds to the replacements of Hx — Hy 



mean freeze-out time tq. The parameters ax, ay and Ht, ax 



n, r. 



y 



and RxRv 



R, 



The azimuthal angle (0) dependence of the invariant 
momentum distribution, eq. (jlOp . can be re-expressed us- 
ing a Fourier-expansion in 0: 



Ni{p) = N^ipup.) 



1 + 2 ^ Vn cos{nq 



(13) 



where Vn are the flow coefficients, in particular V2 is the 
elliptic flow. Sine terms do not appear in the expression 
due to mirror symmetry with respect to the reaction plane. 
The azimuthally averaged transverse momentum distribu- 
tion Ni{pt,Pz) from Ni{p) of eq. PH]) is 



(14) 
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and the elliptic flow: 



V2{Pt,Pz) 



dcl>Ni{p) cos(20) 



(15) 



At mid-rapidity (pz = 0), in reactions of equal- mass nu- 
clei the terms of odd coefficients in eq. (IT3|) disappear due 
to the symmetry — > 4> + n. In this particular case, this 
observable can easily be expressed analytically if we as- 
sume that V2 ^ Vn if n > A, which we may conclude from 
data [25 where V4 ^ V2 ^ 0.01 or smaller, hence it can be 
neglected. Then the formula in eq. (fT3|) simplifies to 



7Vi(p„py,p, =0) =iVi(pt)[l + 2w2C0s(2(/.)], (16) 



where Pt = pI + and Px = Pt cos 4),py = pt sin 0. If we 

evaluate this equation at two appropriate angles (</> — 0° 
and 4> = 45°) then the transverse momentum dependence 
of the coefficient of the elliptic flow can be re-expressed as 
follows 



where we have introduced Teg, the effective slope of the 
azimuthally averaged single particle pt spectra as the har- 
monic mean of the slope parameters in the in-plane and 
in the out-of-plane transverse directions. 



1 _ 1 / 1 ^ 1 



(23) 



The result for the elliptic flow is the following simple 
scaling law: 



V2 



(24) 



where In{z) stands for the modified Bessel function of the 
second kind, /„(z) = (I/tt) exp{z cos9) cos{n9)d9. The 
scaling variable 



w 



pI_ 

4:1711 \Ty T,r. 



(25) 



V2{pt) = ^ 



A^i(Pi,0,0) 
/v, f 0^ 



(17) 



This leads to a simple analytic derivation of V2{pt) at mid- 
rapidity. Numerical investigations at the physical values 
of the model parameters indicate that the formulae have 
about 1% relative error, only. Note that the azimuthally 
averaged transverse momentum distribution in this case 
takes also a simple form as 



Niipt)=Ni 



Pt Pt 







(18) 



However, we can use another but less precise method, 
as well, to find analytic approximation for the elliptic flow. 
In the next procedure a scaling variable w is introduced. 
We show that after some approximations, V2 depends on 
any variable through this variable w only, hence V2{w) is 
a universal function as already pointed out in refs. [4l[T0l 
[^[5] . Both methods were tested against data but the pre- 
vious one proved to describe them with better confidence. 

If we evaluate B(x,p) at mid-rapidity in the limit, 
where the saddle-point coordinates are all small, we get: 



In B{x,p) 



9 

pi 



Pv 



'^-S^, (19) 



2mtTx 2mtTy 2mtTo Tq Tq 
where the direction dependent slope parameters are 



Tx =To + mt R 
Ty = To + mt R. 



b2 



To 



Tq + mtal ■ 
2 ^0 



(20) 
(21) 



The result for the (azimuthally integrated) transverse mo- 
mentum spectrum is: 



Niipt 



9 



{2^y 



\EVC\ . ^ exp 



9 
Pt 



2mtT^ 



cff 



(22) 



This can also be written as 



w = Ek ■ 



T 



cff 



(26) 



where Ek is a relativistic generalization of the transverse 
kinetic energy, defined as 



Ek = 



9 

Pt 

2mt 



and momentum space eccentricity parameter 

T — T 



Tx+Ty 



(27) 



(28) 



In order to compare more easily the Buda-Lund model 
results for the elliptic flow with the azimuthally sensitive 
extension of the Blast Wave model [2^17] , we introduce po 
and p2 so that 



therefore 



Rx = Po(l + P2) 
Ry = Po(l - P2) 

Po = 2 {j^x + 

Rx ~ Ry 

P2 



Rx + Ry 



(29a) 
(29b) 



(30a) 
(30b) 



The elliptic flow in the Buda-Lund model also depends 
on the transverse temperature gradients ax and ay. The 
difference between them actually moderates the difference 
between Tx and Ty and modifles the elliptic flow at given 
Po and p2- On the other hand, at larger momenta and 
small temperature gradient, the following approximative 
proportionality holds: 



V2 OC 



P2 

{Pl-l? 



(31) 
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In the following we detail the results on HBT radii. We 
assume that the fireball is not tilted in the reaction plane, 
i.e. there is no x-z correlation in the emission points of 
the hadron^. There is, however, an angle ip between the 
main axis of the ellipsoidal cross-section of the fireball and 
the outward and sideward axes given by the momentum 
of the hadrons. The former are determined by the orien- 
tation of the reaction plane while the latter are defined so 
that the outward axis agrees with the direction of the av- 
erage transverse momentum of the pair and the sideward 
axis is perpendicular to it. Following eq. (128) in ref. [13 , 
the formula for the two-particle Bose-Einstein correlation 
function can be expressed as 



C2(ki,k2) 



A^2(kl,k2) 



l + Kf2{q»)exp{-qfR^), 



iVi(ki)iVi(k2) 

(32) 

where q = ki - k2 = {q=,q\\,qx,qy), R = [At, R\\, R^, Ry) 
and in the exponential we use the Einstein summation rule 
over the same indices. The longitudinally boost invariant 
temporal^ parallel, sideward, owtward relative momen- 
tum components are defined as 



q= = qo cosh(77) - q^ sinh(77). 




(33a) 


(7|| (72Cosh(?7) - gosinh(ry). 




(33b) 


Qout = {qxKy - qyK.^)/ ^K1 + 




(33c) 


Qside = (qxKa: + qyKy) / ^ 




(33d) 


= -g^^M = + 91 + 'Zx 


+ ql 


(33e) 


= -gi + + ?Le + <?o«t> 




(33f) 



where K — 0.5(fci + k2). The pre-factor Q induces oscil- 
lations within the Gaussian envelope as a function of q\\ . 
This factor is given as 

i7(g||) = cos^ (g|| Zir;) +sin2(q||i?ll Z\r7)tanh^ (77). (34) 

However, it can be approximated by a Gaussian, too, hence 
we may merge it with the longitudinal and temporal emis- 
sion widths as: 



i?,|c = i?|(l + Z\^Vcosh^(^)), (35a) 
At% = smh^{rj)Rl(. + coah^ At'^ . (35b) 



A widely used parameterization of the correlation func- 
tion is given in the longitudinally co-moving frame (LCMS 
[H] , /?/ = 0) of the particle pair and in the out-side- long 
system p5] of Bertsch-Pratt (where = {K^, Kout,0, 0)): 



C2(ki, k2) ~ 1 + A* eyiY>{-qtqjR%), 



(36) 



where i, j — {out, side, long) — (o, s, I). Then the formulae 
for these HBT radii can be expressed by the following 



^ This is justified at RHIC and higher energies, -& — 0. A tilt 
with i9 > may appear at lower energies. 



transformations 
— 2 



R„ 



Ri 



Rl cos^ (p + rI sin^ If + PIAt'^ 



Rx + Ry 



G 

Ry — Rj. 

—5^- cos(2(^) 



(37a) 



—2 



Rj. sin if + Ry cos ip 

Rx+ Ry , - Rx 



cos(2(/j) 



— 2 — 2 

Ry~Rx 



Rf 



cosh^ (77)^1 P + sinh^(?7)Zir^ 
2 

— /3o sinh(^) cosh(^)(i?|| 



= 0. 



(37b) 



(37c) 
(37d) 
(37e) 
(37f) 



4 Comparison to experimental data 

We have determined the best parameter values by fitting 
the analytic expressions for the observables, given in the 
previous section, to experimental data, with the help of 
the CERN Minuit fitting package. Data from 20-30% cen- 
trality class of 200 AGeV Au-|-Au collisions provided by 
PHENIX .26.,27) and STAR [^[^ were used in the analy- 
sis. The fits were performed simultaneously to azimuthally 
integrated transverse mass spectra of positive and nega- 
tive pions, kaons, and (anti)protons |26| , the transverse 
momentum dependence of the elliptic flow parameter V2 
of pions [28^ and to the HBT radii due to pion correla- 
tions as functions of transverse mass and the azimuthal 
angle The results are plotted in Figs. 1-5. 

The interpretation of the model parameters is summa- 
rized in Table [TJ However, for a better understanding the 
results three alternative parameters are introduced that 
can be expressed from the previous ones the following way: 



^sx 



Ri 



2R„ 



2Rl 



T, = 



y 

To 



1 



(38a) 
(38b) 
(38c) 



The two radii parameters correspond to the thermal sur- 
face sizes where the temperature drops to To/2, and the 
parameter Tg corresponds to the temperature of the cen- 
ter after most of the particle emission is over (cooling 
due to evaporation an expansion). Sudden emission cor- 
responds to Te = To, and the At — ^ limit. Also note 
that we use fiB, baryochemical potential. This is calcu- 
lated from the chemical potential of protons and antipro- 
tons: iJ,B = 1/2 (/io,p ~ A*o,p)j as written in Table [TJ 

In Table [21 we present the model parameters obtained 
from simultaneous fits to the data sets. For comparison, re- 
sults are shown from our earlier analysis of 0-30% central- 
ity collisions [30) . too, that was performed with a previous 
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Fig. 1. Buda-Lund model fits to RIfIC 200GeV Au+Au data 
of azimuthally integrated transverse momentum spectra of neg- 
atively charged particles data [26] . 
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Fig. 3. Buda-Lund model fit to RHIC 200GeV Au-f Au data 
on V2 elliptic flow of pions [28] . 
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Fig. 2. Buda-Lund model fits to RHIC 200GeV Au-|-Au data 
of azimuthally integrated transverse momentum spectra of pos- 
itively charged particles [55]. 
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Fig. 4. Buda-Lund model fits to RHIC 200GeV Au-fAu 
data '29' of HBT radii as a function of transverse mass for 
different azimuthal angles. 



version of the model corresponding to the axially symmet- 
ric limit of the current ellipsoidally generalized Buda-Lund 
hydrodynamic model. 

The general observation is that the Buda-Lund model 
parameters describing the source of non-central reactions 
are usually slightly smaller then those of more central colli- 
sions. However, the changes are usually within 2 standard 
deviations, therefore the above statement is based on the 
tendency of the parameters, and on some lower energy re- 
sults not shown here but presented in ref. [30], too. For 
example, the central temperature in these particular non- 
central reactions is below that of the more central ones. 
Also, the transverse geometrical radii at the mean emis- 
sion time are considerably smaller compared to the more 
central values. Moreover, the geometric shape evolution 
due to the asymmetric particle transverse flow in-plane 
(x) an out-of-plane (y) directions results in a source more 
elongated in-plane. Due to the smaller longitudinal source 
size, the parameter corresponding to the formation of hy- 
drodynamic phase is about 10% smaller than that in more 
central colhsions, ro(20 - 30%) = 5.4 ± 0.1 fm/c. The 
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Fig. 5. Buda-Lund model fits to RHIC 200GeV Au-|-Au 
data [29] of HBT radii as functions of azimuthal angles for 
different values of transverse mass. 
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Buda-Lund 


parameter description 


JJcLi cLllic LfcJi b 


To 


Temperature in the center 


a-T- 


Temperature gradient in proper time r 


ax 


Temperature gradient in direction x 




Temperature gradient in direction y 


At 


Distribution width in proper-time r 


Ari 


Distribution width in space-time rapidity r; 


Rx 


Geometrical size in direction x 


Ry 


Geometrical size in direction y 


Hx 


Expansion strength in direction x 


Hy 


Expansion strength in direction y 




Mean freeze-out proper-time 


/is 


Baryochemical potential, 1/2 (^io,p ~ Mo,p) 



Table 1. Description of the parameters of the model. The 
values are meant at the mean freeze-out propoer-time tq. 



elongation in longitudinal direction is similarly smaller, 
Arj{2Q — 30%) = 2.5 ± 0.3. In both cases the baryochemi- 
cal potential is found to be at a low level with respect to 
the proton mass. We emphasize again that the observa- 
tions are based on all the fit results in ref. |30| . 

Note that some of the azimuthally sensitive data have 
large systematic errors that affect the success of fits which 
we had to take into account. The reason for that is the 
difficulty of precise determination of the event reaction 
plane the data are relative to. Several methods are used 
by the experiments to overcome it and we mention those 
applied for the selected data. 

The data set we used for fitting V2 was calculated by 
the four-particle cumulants reaction plane determination 
method that is based on calculations of iV-particle corre- 
lations and non-flow effects subtracted to first order when 
N is greater than 2. The higher N is the more precise 



Buda-Lund 
parameters 


Au+Au 200 GeV 
central (0-30%) 


Au-f Au 200 GeV 
non-central (20-30%) 


To [MeV] 
Te [MeV] 
^iB [MeV] 


196 ± 13 
117 ± 11 
31 ± 28 


174 ± 6 
130 ± 6 
27 ± 16 


Rx [fm] 
Ry [fm] 
Rsx [fm] 
Rsy [fm] 
Hx 

Hy 


13.5 ± 1.7 

Rx 

12.4 ± 1.6 

Rsx 

0.119 ± 0.020 

Hx 


9.5 ± 0.5 
7.0 ± 0.2 

12.8 ± 0.8 

16.9 ± 1.6 
0.158 ± 0.002 
0.118 ± 0.002 


TQ [fm/c] 
At [fm/c] 
Ar] 


5.8 ± 0.3 
0.9 ± 1.2 
3.1 ± 0.1 


5.4 ± 0.1 

2.5 ± 0.2 
2.5 ± 0.3 


XVNDF 


114/208 = 0.55 


269.4/152 = 1.77 



Table 2. Source parameters from simultaneous fits to 
PHENIX and STAR data of Au-I-Au collisions at = 200 

GeV, as given in Fig. 1-5, obtained by the Buda-Lund model. 
For non-central data the value of x^/NDF refers to fits with 
statistical errors, only. 



the event plane determination is, as expected. STAR pub- 
lished two-particle cumulants V2 data in the same refer- 
ence, too, but because of the visible deviations between 
the two kinds of data sets and with respect to the com- 
ments above we used V2{'i} data, only. For further details, 
see ref. ^5] . 

In case of azimuthally sensitive correlation radii, STAR 
has cast about 10% possible systematic errors on the data 
on average. The most likely deviations were assumed to 
take effect on the 'side' and 'out' radii of transverse mo- 
mentum of 0.2 GeV/c. The xVNDF for the fuU fit, includ- 
ing HBT radii with their statistical errors is 269.6/152, 
which corresponds to a very low confidence level. But, 
when we tested our fits with the above mentioned two 
radii of 'side' and 'out' of transverse momentum of 0.2 
GeV/c shifting them within their systematic errors (about 
± 5%) we could achieve an acceptable 1% confidence level 
for the full simultaneous fit. However, without the contri- 
bution of the HBT radii to x^/NDF the confidence level 
is of an acceptable level of 5.1%. 

For comparison to results in refs. [22117] of the az- 
imuthally sensitive extension of the Blast Wave model, 
we have calculated the derivative values given in eqs. (I30al 
[30b)) as po = 1.16 ± 0.05 and p2 = 0.29 ± 0.02. 

5 Conclusions 

We have presented the extension of the Buda-Lund hy- 
drodynamic model from central ultra-relativistic heavy 
ion collisions to peripheral ones. Spectra and the ellip- 
tic flow of identified particles have been described along 
with the azimuthal dependence of two-particle correlation 
function radii in the ellipsoidally symmetric generalization 
of the model. Theoretical predictions were tested against 
RHIC data. From model fits to data of 20-30% central- 
ity class at mid-rapidity source parameters characteriz- 
ing these non-central ultra-relativistic heavy ion reactions 
were extracted. 

The results of our analysis indicate that the central 
temperature in 20-30% centrality reactions is lower than 
that in more central collisions, Tq = 17A ± 6 (stat) MeV. 
In an earlier analysis we showed that in more central (0- 
30%) reactions of the same collision energy this value was 
To = 196 ± 13 (stat) MeV. We have found that the trans- 
verse flow is stronger in the reaction plane than out of 
plane with Hubble constants = 0.158±0.002 and Hy — 
0.118 ± 0.002. The almond shape of the reaction zone ini- 
tially elongated out of plane gets slightly elongated in the 
direction of the impact parameter by the time the particle 
emission rate reaches its maximum. The effect is reflected 
by the geometrical radii in the two perpendicular direc- 
tions at that time, i?2;(in — plane) = 9.5 ± 0.5 (stat) fm, 
i?;^ (out - plane) = 7.0 ± 0.2 (stat) fm. This is the first 
time that an in-plane extended source has been recon- 
structed from simultaneous hydrodynamic model fits to 
identified particle spectra, elliptic flow and azimuthally 
sensitive HBT data in 200 AGeV Au+Au collisions at 
RHIC. 
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The qualitative agreement between model and data 
is apparently good as can be judged from the figures. Al- 
though, in these fits data were used with statistical errors, 
only, which resulted in x^/NDF= 1.77. 
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